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ABSTRACT 



We study four young Pulsar Wind Nebulae (PWNe) detected in TeV 7-rays, 
G21.5-0.9, G54. 1+0.3, Kes 75, and GO.9+0.1, using the spectral evolution model 
developed and applied to the Crab Nebula in our previous work. We model 
the evolution of magnetic field and particle distribution function inside a uni- 
formly expanding PWN considering a time-dependent injection from the pulsar 
and radiative and adiabatic losses. Considering uncertainties in the interstellar 
radiation field (ISRF) and their distance, we study two cases for each PWN. 
Because TeV PWNe have a large TeV 7-rays to X-rays flux ratio, the magnetic 
energy of the PWNe accounts for only a small fraction of the total energy injected 
(typically a few xlO"^). The 7-ray emission is dominated by inverse Compton 
scattering off the infrared photons of the ISRF. A broken power-law distribu- 
tion function for the injected particles reproduces the observed spectrum well, 
except for GO. 9-1-0.1. For GO. 9-1-0.1, we do not need a low energy counterpart 
because adiabatic losses alone are enough to reproduce the radio observations. 
High energy power-law indices at injection are similar (2.5 - 2.6), while low en- 
ergy power-law indices range from 1.0 to 1.6. The lower limit of the particle 
injection rate indicates that the pair multiplicity is larger than 10^. The corre- 
sponding upper limit of the bulk Lorentz factor of the pulsar winds is close to the 
break energy of the broken power-law injection, except for Kes 75. The initial 
rotational energy and the magnetic energy of the pulsars seem anticorrelated, 
although the statistics are poor. 

Subject headings: ISM: individual objects (G21. 5-0.9, G54.1-F0.3, Kes 75, and 
GO. 9-1-0.1) — pulsars: general — radiation mechanisms: non-thermal 
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INTRODUCTION 



A pulsar wind nebula (PWN) is create d by the interaction between the pulsar wine 



and the surrounding supernova (SN) ejecta (IRees &: Gunn 



1974 



Kennel fc Coroniti 



mi) 



The PWN is a shocked pulsar wind composed of a relativistic non-thermal electron-positron 
plasma and magnetic fields. As a result, a PWN shines from radio through TeV 7-rays via 
synchrotron radiation and inverse Compton scattering. The observed spectrum of the PWN 
gives us import ant information on the pair production multiplicity k, of the central pulsar 



( Ide Jager 



20071 ). the m agnetization parameter a in the pulsar wind immediately upstream 



the termination shock flKennel &: Coroniti 



the pulsar wind termination shock (e.g. 



1984). and the part icle acceleration process at 



Fang Sz Zhang 



20101 ). It should be noted that 



the radiation spectrum depends on time through the evolution o: 



expansion of ; 



Gelfand et al. 



le PWN, and radiative and a diabatic losses (e.g. 



the energy iniection, the 



Tanaka fc Takahara 



2010; 



2009 



Bucciantini et al. 



20101 ). 



Through the recent development of 7-ray observations, many PWNe have been 



discovered in 7-rays and many detailed struct ures of them have also 



Deen found in other 



wavelengths (e.g., HESS J1640-465 studied by lAharonian et al. 



2006 



Lemiere et al. 



20091) 



Although the observed PWNe show common characteristics, such as flat radio spectrum 



Gaensler fc Slane 



20061 ) 



and X-ray spectral steepening with distance from the pulsar (c.f., 
the differ. We discuss these individual differences of PWNe in this paper. They include the 
age and the spin-down power of the central pulsar, the observed size of the PWN and the 
TeV 7-rays to X-rays flux ratio. Here we study these characteristics because they constrain 
the physical condition of the pulsar magnetosphere and the pulsar wind. 



We built a spectral evolution model of PWNe in our previous work (ITanaka fc Takahara 



2OIOI ). The application of the model to the Crab Nebula well reproduces the current 



observed spectrum and the radio flux evolution. In this paper, we apply this model to four 
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young PWNe detected in TeV 7-rays, G21.5-0.9, G54.1+0.3, Kes 75, and GO.9+0.1. We 
pick up these four PWNe according to three criteria. (1) They reveal observed non-thermal 
spectrum at least in radio, X-rays, and TeV 7-rays. (2) They have a central pulsar with 
known period and its derivative. (3) They have an almost spherical shape with a known 
angular extent, i.e., they are young enough not to reveal signatures of interaction with SNR 
reverse shock. For the application to the PWNe other than the Crab Nebula, we include 
the interstellar radiation field (ISRF) in infrared and optical bands as the target photons of 
the inverse Compton scattering. 

We assume a constant expansion v elocity of the PWN. This assumption is appropriate 



to young PWNe. 



Gelfand et al. 



(j2009f ) considered the dynamical evolution of a PWN 
inside a supernova remnant (SNR). Their model of the PWN expansion gives a more 
realistic dynamical evolution than the constant expansion velocity and also applies to old 
PWNe. However, their expansion model includes many unknown parameters to determine 
the dynamical evolution including the energy of the SN explosion, the mass of the SN 
ejecta and the density of the surrounding interstellar medium. We factor these uncertain 
quantities into one parameter: the constant expansion velocity. Note that a simple 
estimate of the radius of th e PWN in the early phase of its evolution, -RpwN oc t^^^ (e.g.. 



van der Swaluw et al. 



20011 ). is close to an expansion at constant velocity. 



In Section |2l we describe our model and slightly improve it for the application to young 
TeV PWNe other than the Crab Nebula. In Sections [3] -[6l we apply the model to G21.5-0.9, 
G54.1+0.3, Kes 75, and GO.9+0.1 and discuss about their individual characteristics. A 
comparative discussions about the young TeV PWNe we studied, including the Crab Nebula 
and the conclusions are made in Section [71 
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THE MODEL 



We give here a brief 



(ITanaka Sz Takahara 



description of our spectral evolution model of PWNe 



20101) • For application to young TeV PWNe other than the Crab 
Nebula, we improve the target photon fields for inverse Compton scattering. We include 
the ISRF in infrared and optical bands, which were ignored in the previous work. Lastly, 



we describe the fitting procedure. 



2.1. Basic Ingredients of Our Model 

We use a one-zone model of PWNe. A PWN is a uniform sphere expanding at a 
constant velocity fpwN, i-e., the radius of PWN is given by -Rpwn(^) = "^^pwn^- The contents 
of the PWN are the magnetic field and the relativistic non-thermal electron-positron 
plasma and they are injected from the central pulsar. The magnetic energy injection 
Emagit) = rjL{t) and the particle energy injection i?part(^) = (1 — v)L(t) are characterised 
by the fraction parameter r] {0 < rj < 1) and the spin-down power L{t). The evolution of 
the spin-down power L{t) is given by 

n + l 

L{t) = Lo f 1 + -J , (1) 

where Lq is the initial spin-down power and tq is the spin-down time. For t > tq, the 
total energy injected into the PWN up to time t (i?tot(^) = Jq L{t')dt') approximately 
corresponds to the initial rotational energy of the central pulsar Lq-tq = IVL^/{n — 1), where 
VLq = 2tt/Po is the initial angular velocity. For t < tq, E'tot is smaller than Lq ■ tq. We need 
four quantities, the current pulsar period P, its time derivative P, braking index n and the 
age of pulsar tage, to fix the evolution of the spin-down power, assuming that the moment 
of inertia of the pulsar is lO^^g ■ cm^. Note that three time scales appear in the spin-down 
evolution of the pulsar: the age of pulsar tage, the spin-down time tq and the characteristic 
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age Tc- They are related to each other through the simple relation 



TL — 1 

T"c = (^0 + W)- (2) 



The customarily used braking index n = 3 gives Tj. = tq + tage- 

We assume that the distribution of the particles in the PWN is isotropic, and then the 
evolution of the particle distribution A^(7, t) is given by the continuity equation in energy 
space, 

^iV(7,t) + |- (7(7,t)iV(7,t)) = ginj(7,t). (3) 
We assume that the particle injection Qin]{lit) follows the broken power-law distribution 

^ , I Qo(t)(7/7b)-^^ for7min<7<7b, 

<9mj(7,i)=< (4) 

I Qo(t)(7/7b)-^^ for7b<7<7 max ; 

where 7 is the Lorentz factor of the relativistic electrons and positrons. We introduce the 
parameters of the injection spectrum, 7min5 7b5 7max5 Pi < 2 and p2 > 2, which are the 
minimum, break and maximum Lorentz factors and the power-law indices at the low and 
high energy ranges of the injection spectra, respectively. We require that the normalization 
Qolt) satisfies (1 — ri)Lit) = J^™'"' Qinj(7, t)7meC^rf7, where rrie and c are the mass of an 
electron (or positron) and the speed of light, respectively. We consider the cooling effects 
of the relativistic particles 7(7, t) including the synchrotron radiation 7syn(7;^); the inverse 
Compton scattering off the ISRF 710(7) and the adiabatic expansion 7ad(75^)- The detailed 
description of the ISRF will be discussed in Section 12. 2[ The cooling time of the particles 
Tcooi(7,^) = 7/ 17(751^)1 is an important time-scale in addition to tagc and tq, and these three 
time-scales characterize the evolution of the particle distribution. 

For the magnetic field evolution, we assume the following form of the magnetic field 
energy conservation, 

^{RpwAt)f ■ ^-^P^ = [\mdt' = r^E,,,{t). (5) 
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The magnetic field approximately evolves as B{t) cc t ^ for t < tq and B{t) oc t for 



t > tq. Note that some justifications o: 
the previous paper (see Section 2.2 of 



this magnetic field evolut ion model are discussed in 



Tanaka &: Takahara 



2010h 



For the calculation of the radiation spectrum, we assume that the radiation is isotropic. 
The radiation processes which we consider are synchrotron radiation and inverse Compton 
scattering off the synchrotron radiation (SSC) and the ISRF including the cosmic microwave 
background (IC/CMB) radiation, infrared photons from dust grains (IC/IR) and optical 
photons from starlight (IC/OPT). 

Lastly, we mention the evolution of the synchrotron cooling break frequency for 
convenience in the later discussion. Assuming that the inver se Compton cooling is 



ineffective in most of the evolutionary phase (e.g.. Figure 4 of iTanaka fc Takaharal 120101 ). 
we divide the particle distribution into two populations; one corresponds to the high energy 
particles which lose energy by synchrotron cooling, the other is the low energy particles 
which lose energy by adiabatic cooling. The critical energy 7c (t) is defined by equating the 
synchrotron cooling time Tsyn(7,t) with the adiabatic cooling time Tad(7,t) = t, which is: 



For 7 > 7c(t), the synchrotron cooling dominates over the adiabatic cooling and vice versa 
for 7 < 7c (t). The characteristic frequency of synchrotron radiation z/gyn is a function of the 
particle energy and the magnetic field. Because of the rapid decrease of the magnetic field 
strength, the synchrotron cooling break frequency udt) = i^synil c{t), B{t)) increases with 
time, which is: 



- 8 - 



2.2. Galactic Interstellar Radiation Field as the Target of Inverse Compton 

Scattering 



We use the ISRF model described by iPorter et al.l (120061 ) as a reference. Their ISRF 



model depends on the distance from the galactic center r and the height from the galactic 
plane z. Although the contribution of localized dusts and stars to the ISRF changes the 
actual ISRF around the objects, we ignore the possibility of local effects, except for the 
model 2 of G54.1+0.3. 

We assume that the ISRF has three components. The CMB is a blackbody radiation 
with a temperature Tcmb = 2.7K. The spectra in infrared and optical bands are 
modified blackbodies characterized by a temperatures and an energy density (Tjr, Uyr) 
and (TopT, t^oPT), respectively. Because the temperature dependence of the IC/IR 
(TiR ~ 30 - 50ir) and the IC/OPT (Topt ~ 3500 - 4000/s:) is weak, we fix Tir = 40K 
(mildly Klein-Nishina regime), and Tqpt = 4000K (mostly Klein-Nishina regime) for all 
objects. On the other hand, the IC/IR and the IC/OPT luminosities strongly depend on 
the energy densities f/jR and t/oPT which vary with galactic locations. From the distance 
to each PWN and its galactic coordinates, we can roughly estimate (r, z) for each object. 
Adopted values for each object are listed in Table [H 



2.3. Fitting Procedure 

Here, we summarize the model parameters and describe the fitting procedure. We 
adopt the distance to the object c?, the angular extent of the object, and the central pulsar 
parameters P, P and n from observations. Once the distance to the object is fixed, we can 
determine the radius of the object -RpwN and the energy densities of the ISRF f/jR and 
f/oPT- 
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Fitting parameters are the age tage (assuming that the age of the central pulsar is the 
same as that of the PWN), the fraction parameter rj, and the parameters of the injection 
spectrum in Equation (jl]), 7max5 Ih, Imin, Pi and p2- In our spectral evolution model, the age 
tage and the fraction parameter 77 of PWNe are primarily determined by the absolute power 
of the synchrotron radiation and the po' wer ratio of the inverse Com pton scattering to the 



synchrotron radiation (see Section 3.3 of iTanaka fc Takaharall2010l ). The parameters of the 



injection spectrum in Equation (jl]) are determined by comparing the calculated spectrum 



with the detailed shape of the observed spectrum (see Section 3.4 of 



Tanaka fc Takahara 



20101 ). The parameters 7min and 7max are relatively unconstrained. We obtain only upper 
limit of 7min from the lowest frequency in radio observations. We choose lower limit of 
7max so that i^syn(7max, B) bccomcs an order of magnitude larger than the observed highest 
frequency in X-rays, since no clear spectral rollover has been observed for the four PWNe. 

The fitted parameters which reproduce the multi-wavelength observations are almost 
uniquely determined except for 7inin and 7max, if we fix the distance to the PWNe and the 
energy densities of the ISRF. To understand effects of uncertainties in the observed flux, 
measured distance and the ISRF model on our calculations, we study two different models 
for each PWN in Sections [3] - [61 

The adopted and fitted parameters gives the values of the expansion velocity fpwN, the 
current magnetic field strength of the PWNe -Bnow and the parameters of the central pulsars 
To and Lq ■ tq. The adopted, fitted and derived values for each model are listed in Tabled! 



G21. 5-0.9 



G21.5-0.9 is a com posite SN R and its PWN 



1989 



Bock et al 



2OO1I), infrared flGallant fc Tuffs 



is ob served in radio (ISalter et al 



19981 ). X-rays (ITsujimoto et al. 



2010; 
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De Rosa et al. 



20091 ) ■ and TeV 7-rays ( iDjannati-Atai et al. 



2003). In X-rays, G21.5-0.9 



consists o f an outer halo ^ 150^^ in rad i us, which is the SNR component, and a PW^ » 



in radius (IMatheson fc Safi-Harb 



20051: 



Bocchino et al 



2005 



Matheson &: Safi-Harb 



40" 



2010h . 



In radio, the outer halo is not o 



in the X-rays (iBock et al 



2001 



jserved and the size of the PWN is comparable w ith that 



Bietenholz &: Bartel 



2008 



Bietenholz et al. 



2010). Infrared 



observations of PWNe are difficult and G21.5-0.9 is one of a few PWNe whose non-thermal 



infrared spectrum is observed (iGallant &: Tuffs 



the flux of the PWNe component in GeV 7-rays (lAckermann et al 



19981). Fermi LAT puts an upper limit on 



20111 ). The central 



pulsar of G21.5-0.9 (PSR J1833-1034) is observed in radio and GeV 7-rays with a period 
P = 6.19 X lO^^sec, its ti me derivative P = 2 .02 x 10~^^sec ■ se c~^ and unknown braking 



index (we assume n = 3) (ICamilo et al. 



of the central pulsar is 4.9kyr, but 



2006 



Abdoeta 



20091) . The characteristic age Tc 



Bietenholz fc Bartell ( 120081 ) suggeste d t^g^p ^ 900yr frorn 



the observation of the expansion rate of the PWN. On the other hand, Wang et al. 



suggested that G21.5-0.9 might be the historical sup ernova in 48 BC, i.e ., t. 



age 



(119861) 
2kyr. 



20081 ) and then we 



We assume that the distance to G21.5-0.9 is 4.8kpc ( iTian fc Leahyl 
approximate G21.5-0.9 PWN as a sphere of radius l.Opc at (r, z) ~ (4kpc, 80pc). Lastly, we 
get (t/iR, f/opx) = (l.OeV/cm"^, 2.0eV/cm^) as the energy density of the ISRF. Considering 
that the observed infrared spectrum of G21.5-0.9 is uncertain, we investigate two cases with 
(model 1) and without (model 2) infrared fltting. 



3.1. Model 1 

Figure [1] shows the model spectrum of G21.5-0.9 with the observational data including 
the infrared observation. We flt the data with the parameters 77 = 1.5 x 10~^, tage = l.Okyr, 
7max = 2.0 X 10^ 7b = 1.2 X 10^ 7i„in = 3.0 X 10^ pi = 1.0, and p2 = 2.55. The fltted 
fraction parameter is three times larger than the Crab Nebula, but still much smaller than 
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unity. The fitted age tage = l.Okyr deviates from tlie cliaracteristic age of tlie central pulsar 
Tc = 4.9kyr. Because the ratio of the fitted age to the characteristic age of the pulsar is 
described as tagc/Tc = 2/(r2 — 1) ■ [1 — {Pq/ PY'"'-^^''], the current pulsar period is almost the 
same as the initial one {Po/P ~ 0.9). The derived expansion velocit y fpwN = 980km/sec 



and t he fitted age are consistent with the observed expansion rate ( iBietenholz &: Bartel 
20081 ). From the fitted age, we obtain a spin-down time tq = 3.9kyr and an initial rotational 
energy of the pulsar Lq ■ tq = 6.5 x lO^^erg. The current total energy injected into the 
PWN _Etot(lkyr) = 1.3 x lO'^^erg is significantly below Lq ■ tq, since tq > tage- The current 
magnetic field strength of G21.5-0.9 turns out to be -Bnow = 64/iG. We choose 7min and 7max 
by taking z/syn(7min, ^now) ~ lO^Hz and z/syn(7max, ^now) ~ 3 X 102°Hz, respectively. 

There are two breaks in the calculated synchrotron spectrum in Figure [1] One is the 
synchrotron cooling break frequency ~ 3 x lO^^Hz, which is determined from Equation 
([7]), the other is the break in the particle injection z^syn(7b5 -Bnow) ~ 8 x lO^^Hz, which is 
much smaller than Uc- Since 2 < p2 < 3, the synchrotron spectrum in Figured] is peaked at 
around Uc- Focusing on the calculated TeV spectrum in Figure [H the IC/IR is found to 
dominate and the SSC contribution is negligible. Although the energy density of the ISRF 
in optical band is twice as large as that in infrared band, Klein-Nishina effect significantly 
reduces the IC/OPT. The calculated flux of the GeV 7- ray is almost an order of magnitude 



Ackermann et al 



(1201 ih . In hard 



below the upper limit of the PWN component given by 
X-rays, the observed spectrum in 5 x 10^*^ - 2 x lO^^ Hz is well reproduced , while the 



pulsar emission appears to dominate above 2 x lO^^Hz (iDe Rosa et al. 



2OO9I ). The observed 



spectrum in soft X-rays (< 3 x 10^'^Hz) and the observed sharp break around 3 x lO^^Hz 
are difficult to reproduce in our model. The observed soft X-ray spectrum is much harder 
than the calculated spectrum and such a hard spectrum does not smoothly connect to 
the observed infrared spectrum. We will discuss about these discrepancies between the 
calculated and observed spectra in soft X-rays in Section | 



- 12 - 



Figure [2] shows the evolution of the emission spectrum (left panel) and that of the 
particle distribution (right panel) of G21.5-0.9 with the use of the same parameters as 
in Figure [TJ While the synchrotron flux decreases with time, the inverse Compton flux 
increases. This feature is due to the decrease of the magnetic fleld strength with time 
and to the increase of the particle number as seen in the right panel. Because the energy 
injection from the pulsar will continue till the tim e t ^ Tn = 3.9kyr, the part icles in the 



PWN increase (see the discussion in Section 3.2 of iTanaka &: Takaharal (|2010[ ) in details). 
In the left panel, we can see the evolution of the synchrotron cooling break frequency t'c(^) 
(from ~ 7 X lO^^Hz at 300 yr to ~ lO^^Hz at 10 kyr) as is predicted by Equation (^^. 



3.2. Model 2 

The calculated spectrum in model 1 does not reproduce the observed spectrum in 
infrared and soft X-rays at the same time. Considering the uncertainties in the infrared 
observations, we try to flt the observed spectrum ignoring the infrared band. Figure [3] 
shows the model spectrum of G21.5-0.9 without infrared fltting. All the fltted parameters 
are similar to model 1 (see Tabled]), but a slightly smaller value of the fraction parameter 
7] = 8 X 10^^ (Aiow = 47/uG) is allowed. Consequently, the synchrotron flux in infrared is 
smaller than model 1 and the synchrotron cooling break frequency z/c ~ 7 x lO^^Hz is a 
little larger than model 1. The observed spectrum in soft X-rays and the observed sharp 
break at around 3 x lO^^Hz are not well reproduced, even if we ignore infrared observations. 
The hard X-ray observation (5 x 10^^ - 2 x lO^^Hz) is reproduced as well as model 1. 
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3.3. Discussion 

Both model 1 and 2 reproduce the observational data reasonably well. The fraction 
parameter rj, i.e., the magnetic field strength B^^w, is the only notable difference between 
the models. Future observations in infrared and optical bands will provide more accurate 
spectrum of G21.5-0.9 and could distinguish the models. As for the other parameters, the 
values of 7b and pi are smaller than the Crab Nebula in both models. However, the high 
energy power-law index at injection p2 is very similar to the Crab Nebula. 



de Jager et al 



(120091). They solved 



We compare our results with the model by 
hydrodynamic equations and the evolution of the magnetic field is separately calculated 
by the induction equation. The evolution of the broad band spectrum is calculated with 
a one-zone approach with the use of the volume averaged magnetic field. Their adopted 
energy density of the ISRF in infrared band f/jR = l.OeV/cm^ and tage = Ikyr are the 
same as ours. They obtain a particle injection distribution with pi = 1.0, p2 = 2.6 
and 7b ~ 8 X 10^, which are almost the same as ours. Their obtained current magnetic 
field strength ~ 24/iG is smaller than ours, -Bnow = 64yuG and 47/iG for model 1 and 2, 
respectively. Their results are basically consistent with ours. 

Lastly, we discuss the di s crepa ncy between our model spectra and the observations in 
soft X-rays. jTsujimoto et al.i (|201Q) show that there is an observed spectral break between 
the soft and hard X-rays at ~ 3 x lO^^Hz. One might consider that this spectral break 
corresponds to the synchrotron cooling break frequency Uc. With t^ge = Ikyr, the value 
of z/c ~ 3 X lO^^Hz requires the magnetic field strength to be ~ 6/iG from Equation ([7]). 
However, the current magnetic field of ~ 6/iG seems unlikely because the observed 7-rays 
to X-rays fiux ratio would demands the local ISRF energy density around G21.5-0.9 to be 
much smaller than the contribution from the CMB. Note that it is hard to reproduce the 
observed sharp spectral break at ~ 3 x lO^^Hz, even if the magnetic field strength is ~ 6yuG. 
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This is because the evolution of the magnetic field B{t) and the particle injection Qmj(7,^) 
makes the synchrotron cooling break rather smooth. We thus conclude that any one- zone 
spectral evolution models fail to reproduce the current soft X-ray observation. One possible 
resolution is a modification of the observed spectrum in soft X-rays with relat ively large 



interstellar extinction toward G21.5-0.9 as discussed by 



Safi-Harb et al.l (120011 ). Another 



po ssibility is t 



by 



l e effec t of spatial variation of the soft X-ray photon index Fx as discussed 



Slane et al. 



(120001 ). 



4. G54.1+0.3 



G54.1+0.3 is a center filled or possibly a composite SNR and its PWN component is 



obser ved in radio ( 



2010h . X-rays 



Green 



Lu et al 



1985 



Velusamy et al 



1986 



Velusamy fc Becker 



1988 



Lang et al 



200 ll ). and TeV 7-rays (iAcciari et a. 



jet-torus structure is observed similarly to the Crab Nebula (iLu et aL 



.il2010f). In X -rays, a 



of the radio PWN ^ 2.5' x 2.0' is comparable with that in X-rays (ILang et al. 



2002). The size 



Bocchino et al. 



by 



Bocchino et al 



2OIOI : 



20101). R ecently, diffuse X-ray emission surrounding the PWN was detected 



2OIOI ). which is a po ssible counterpart of SN ejecta component, while a 



possible SNR shell is observed in radio ( 



has been observed in infrared (iKoo et al. 



ang et al 



2008 



2010). T 



le re gion around G54.1+0.3 



Temim et al.l 



2010l ). Although no infrared 



counterpart of G54. 1+0.3 PWN is observed, there are several bright sources around 
G54. 1+0.3. The central pulsar of G54.1+0.3 (PSR J1930+1852) is observed in radio and 
X-rays with a period P = 1.36 x lO^^sec, its time derivative P = 7.51 x lO^^^sec ■ sec~^ 



Tc = 2.9kyr) and unknown braking index (jCamilo et al.ll2002l ). So that we assume n = 3. 



he actual age of G54.1+0.3 is unknown. We assume a distance to G54. 1+0.3 of 6.2kpc 



(ILeahy et al. 



20081 ) and then we assume that G54. 1+0.3 PWN is a sphere of radius 1.8pc at 
r, z) ~ (7.5kpc, Opc). We get (f/m, f/opx) = (0.5eV/cm^, 0.5eV/cm^) for the energy density 
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of the ISRF in model 1. We also investigate in model 2 the probability that the observed 
infrared sources around G54. 1+0.3 significantly contribute to the local ISRF. 

4.1. Model 1 

Figure m shows the model spectrum of G54.1+0.3 for model 1 with the observational 
data. We can reproduce the observed non-thermal spectrum with the parameters 
r/ = 3.0 X 10-^ 4ge = 2.3kyr, 7^,, = 1.0 x 10^, 7b = 3.0 x 10^ 7min = 2.0 x 10^ 
Pi = 1.2, and p2 = 2.55. The fraction parameter is 0.06 times the Crab Nebula's value. 
The fitted age tage = 2.3kyr is comparable with the characteristic age of the central pulsar 
Tc = 2.9kyr, and the current period is almost twice as large as the initial period. The 
corresponding expansion velocity fpwN = 770 km/sec is less than a half of that of the 
Crab Nebula 1800 km/sec. We obtain a spin-down time tq = 0.6kyr < tage and an initial 
rotational energy Lq ■ tq = 5.4 x lO^^erg. The current total energy injected into the PWN 
i?tot(2.3kyr) = 4.3 x lO^^erg is close to the value of Lq ■ tq. The current magnetic field 
strength of G54.1+0.3 turns out to be -Bnow = 6.7/iG, which is much smaller than that of 
the Crab Nebula and G21.5-0.9. 

The synchrotron cooling break frequency z/c ~ 5 x 10 ^^Hz is much larger than the 
characteristic synchrotron frequency corresponding to the break energy z^syn(7b5 -Bnow) ~ 
7 X lO^^Hz. In 7-rays, the ISRF energy density in infrared band is twice as large as that of 
the CMB, but the contributions of the IC/CMB and the IC/IR are comparable because the 
IC/IR is in mildly Klein-Nishina regime. Going into details, the observed X-ray spectrum 
seems consistent with the model spectrum, while the observed 7-ray spectrum seems softer. 
Because particles with almost the same energy contribute both to the observed emission in 
X-rays (SYN) and 7-rays (IC/CMB in Thomson regime), it is difficult to fit the spectral 
slopes in both frequency ranges at the same time. 



- 16 - 



Figure [5] shows the evolution of the emission spectrum (left panel) and the particle 
distribution (right panel) of G54.1+0.3 with the use of the same parameters as in Figure 
m Both the synchrotron flux and the IC flux decrease with time, but the IC flux deceases 
more slowly. This feature is because the decrease of the magnetic field strength is 
much faster than the particle distribution while the ISRF energy density is constant. 
This behavior of the parti cle distribution has already been discussed in Section 3.2 of 



Tanaka &: Takaharal (l2010l ). In short, the evolution of the particle number is approximately 



expressed as A^(7, t) Q{'y,t) ■ t ■ Tcoo\{i, t)/(t + Tcooi(7, t)) and the injection from the pulsar 
decreases as Q{'-f,t) oc after the time t ~ tq = 600yr. In synchrotron cooling regime 
(tcooI = Tsyn < Tad ~ t) , the partlclc number increases as oc t, while in adiabatic cooling 
regime (tcooi = Tad), the particle number decreases as oc t~^, where we use Tgyn oc for 
t > tq. The left panel of Figure [S] shows that the synchrotron cooling break frequency 
z^c(^) will continue to increase with time till an age of 3kyr (Equation ([7j)). Because 
7c(10kyr) > 7max, the synchrotron spectrum will be peaked at z/syn(7max), not Uc, at an age 
of lOkyr. 



4.2. Model 2 



The region ar ound G54. 1+0.3 has been observed in infrared (iKoo et al. 



2008 



Temim et al 



20101 ) and observations suggest that the ISRF around G54. 1+0.3 could be 
larger than that of the average values of the Galaxy adopted in model 1. Figure [6] shows 
the model spectrum of G54.1+0.3 for model 2 when the energy density of the ISRF in 
infrared band has been chosen to be four times larger than model 1 ([/jr = 2.0eV/cm^). 
All the fitted parameters are similar to model 1 except for the fraction parameter and 
the fitted age; they are r] = 2.0 x 10"'^ and tage = l-7kyr, respectively (see Tabled]). 



The fitted fraction parameter is close to the values of the Crab Nebula. The fitted 



-17- 



ar to the i ndepe ndent estimate of 1.5 kyr from the dynamical 



age tage = 1.7kyr is simi 

interaction with SNR by IChevaherl (|2005[ ). Because the fitted age tage changes from model 
1, the derived parameters change also from model 1, wpwN = 1040km/sec, tq = 1.2kyr and 
Lq ■ tq = 2.6 X lO^^erg (i?tot(l-7kyr) = 1.6 x lO^^erg), respectively. The current magnetic 
field strength -Bnow = lOyuG becomes a little larger than model 1. We can see the trend 
that a larger energy density of the ISRF leads to a larger r] and -Bnow and a smaller value 
of -Etot(tage) a nd ta^p. This trend c a n be understood following a similar discussion made in 



Section 3.3 of 



Tanaka fc Takaharal ( |2010| ) (see Section [7?!]) . 



4.3. Discussion 



We favor model 2 because the observed 7-ray spectrum is somewhat better reproduced 
by model 2, although both model 1 and 2 can reproduce the observational data reasonably 
well; model 2 has a softer 7-ray spectrum than model 1 through the contribution of the 
IC/IR is larger than the IC/CMB. Because the difference of the ISRF energy density in the 
models appears at the frequency where the Klein-Nishina effect works, future observations 
in higher energy 7-rays (1 - lOOTeV), such as CTA, would provide better information on 
the correct values of the ISRF. In contrast to 77 and tage, the parameters of the particle 
injection for both models are similar. The low energy power-law index at injection pi is 
different from the Crab Nebula and G21.5-0.9, but the high energy power-law index at 
injection p2 is very similar. 



The spectral evolution of G 54. 1+0.3 was a lso studied by 



Li et al. 



pure-lepton model (Figure 2 of 



Li et al. 



of the broad band spectrum (Figure 3 of 



fcoioh . In addition to 



20101). they stu died a lepton-hadron hybrid model 



Li et al. 



20101). In the pure-lepton model, they 



obtained the parameters pi = 1.2 and 7b = 5 x 10^, which are similar to ours. However, 
their obtained p2 = 2.8 and adopted Ui^ ~ 3.3eV/cm^ are different from ours in model 2, 
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P2 = 2.55 and f/jR ~ 2.0eV/cm^, respectively. Difference appears in the X-ray spectrum, 
which in their model seems softer than ours. The current magnetic field strength in their 
model ~ lOyuG is the same as our model 2. We conclude that their pure-lepton model is 
almost consistent with model 2 of ours except for the X-ray spectrum. 



Li et al. 



( I2OIOI ) argued that the lepton-hadron hybrid model {B ~ 80/iG) is better than 



the pure-lepton model because the current magnetic field B 



is much weaker than the observational indication by 



Lang et al. 



OuG in pure-lepton model 



fl2010h . 



Lang et al. 



fl2010h 



estimated an equipartition magnetic field of 38 from the radio luminosity of PWN and a 



magnetic field of 8 - 200 uG 



these estimates by iLang et al. 



r om th e lifetime of X-ray emitting particles. We consider that 



(120101 ) are not robust. Generally, the magnetic field strength 



is sub-equipartition for all PWNe we stu died, i.e., y « 1, so that B = 10/iG is just a 



reasonably expected value for this PWN. 



Li et al 



(I2OIOI ) also argued that the observed 



7-ray photon index is better fitted by the lepton-hadron hybrid model than the pure-lepton 
model. However, the calculated 7-ray spectrum in the pure-lepton model changes with local 
ISRF energy density and temperature. We believe that leptonic model with small magnetic 
field is consistent with the current observations. 



Kesteven 75 



Kes 75 is a composite SN R and it s PWN component i s observed in radio ( Salter et al. 



1989 



Bock &: Gaensler 



2005), X- 



7-ray s (iDjannati-Atai' et al. 



rays (IHelfand et al. 



2003 



McBride et al. 



2007h , and TeV 



20071 ). In X-rays, a jet-torus structure is observed flNg et al. 



20081 ). This jet-torus structure is surrounded by a diffuse X-ray nebul a, which has almost 



the same extent of the radio PWN ~ 26" x 20^' ( 



SNR shell is observed in radio and X-rays (IBock fc Gaenslei 



Helfand et al. 



2005 



2003h. A part o : 



f the 



Helfand et al 



20031) 



The central pulsar of Kes 75 (PSR J1846-0258) is observed in X-rays, with a period 
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P = 3.26 X 10 ^ sec, its time derivative P = 7.08 x 10 ^^sec ■ sec ^ and braking index 



n 



2.65 ± 0.01 (ILivingstone et al. 



20061 ). Although the characteristic age of the central 
pulsar Tc = 0.7kyr suggests it is very young, the actual age of the Kes 75 is unknown. 
PSR J1846-0258 has a large surface dipole magnetic field and a magnetar-like burst was 



observed in 2006 ( iGavriil et al. 



20081 ). No radi o counterpart of PSR J 1846-0258 is observed, 



common to most of the magnetar cand idates (lArchiba 



cases for the distance to Kes 75: 6kpc ( iLeahy &: Tian 



d et all 120081 ). We consider 



two 



20081) and 10.6kpc fISu et al 



20091 ) 



as model 1 and 2, respectively. We assume that Kes 75 PWN is a sphere of radius 0.29pc 
at (r, z) ~ (4kpc, 30pc) or 0.5pc at (r, z) ~ (5kpc, 50pc) for model 1 and 2, respectively. 
The energy density of the ISRF is (t/m, t/oPT) = (1.2eV/cm^, 2.0eV/cm'^) for model 1 and 
(?7iR, ?7opt) = (l.OeV/cm^, 2.0eV/cm^) for model 2, respectively. 



5.1. Model 1 

Figure [7] shows the model spectrum of Kes 75 when the distance to Kes 75 is 
assumed to be 6kpc together with observational data. We fit the data with the parameters 
= 5.0 X 10-5, = 0.7kyr, 7^ax = 2.0 x 10^, 7b = 2.0 x 10^, 7^;, = 5.0 x 10^, p, = 1.6, 
and P2 = 2.5. The fraction parameter is very small and two orders of magnitude smaller 
than the Crab Nebula. The fitted age tage = 0.7kyr is very close to the characteristic age 
Tc = 0.7kyr and the expansion velocity fpwN = 420km/sec is rather slow. The pulsar 
parameters tq = 0.2kyr and Lq ■ tq = 1.5 x 10^'^erg are also smaller than those of the other 
PWNe. The small tq may be related to the large magnetic field of the PSR J1846-0258. 
The current total energy £'tot(0.7kyr) = 9.1 x 10^''erg is roughly half of Lq ■ Tq. The current 
magnetic field strength of Kes 75 turns out to be -Bnow = 20/iG. Despite r] being more than 
an order of magnitude smaller than in the other PWNe, -Bnow is not so different because the 
size of the Kes 75 PWN is small, and accordingly, fpwN is small. 
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Focusing on the detailed spectral features, the observed flux in hard X-rays is not well 
reproduced and a few times larger than the model prediction although the characteristic 
synchrotron frequency corresponding to 7max extends to ~ 3 x 10^°Hz. The hard X-ray 
observation is difficult to reproduce because the synchrotron cooling break frequency 
z/c ~ 2 X 10^ ''Hz is located at the soft X-rays (see discussion in Section [5731) . The 7-ray 
emission is the IC/IR dominant like G21.5-0.9 and G54.1+0.3. 

Figure [8] shows the evolution of the emission spectrum (left panel) and the particle 
distribution (right panel) of Kes 75 with the use of the same parameters as in Figure [71 
Both the synchrotron flux and the IC flux decrease with time. This feature is similar to 
054.1+0.3 (Figure [5]) because of the small spin-down time tq = 0.2kyr. 

5.2. Model 2 

Figure M shows the model spectrum of Kes 75 when the distance to Kes 75 is assumed 
to be 10.6kpc together with observational data. We calculate the spectrum with the 
parameters i] = 6.0x 10"^, tagc = 0.88kyr, 7jnax = 1.0 x 10^, 7b = 5.0 x 10^, 7min = 5.0 x 10^, 
Pi = 1.4, and p2 = 2.5, but the calculated 7-ray flux is a few times smaller than the observed 
one. The fraction parameter is an order of magnitude smaller than model 1. The fltted age 
tage = O.SSkyr is almost the maximum value given by Equation ([2]). The current magnetic 
fleld strength -Bnow is 24yuG similar to model 1 because the adopted ISRF energy density 
is similar. On the other hand, Etot = 1-7 x 10^'^erg is more than two orders of magnitude 
larger than model 1 given the larger 7-ray luminosity. The pulsar parameters tq = 3yr and 
Lq ■ To = 2.1 X lO^^erg ~ Etot are more extreme than model 1. 

To better reproduce the observed 7-ray flux in model 2, ii^tot larger than 1.7 x 10^'^erg is 
required with the adopted ISRF energy density. i?tot > 1.7 x lO^^erg leads to tq < 3yr and 
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Lq ■ tq > 2.1 X lO^^'erg. However, tq = Syr is extremely short compared with tage = 0.88kyr 
and Lq ■ To = 2.1 X 10^°erg is three times larger than the Crab Pulsar. On the other hand, 
these extreme Tq, Lq ■ tq and also r] may be allowed given the unique properties of the 
central pulsar PSR J1846-0258, and in this case an energy density of the ISRF a few times 
larger would reproduce the observed 7-ray flux. 



5.3. Discussion 

We favor model 1 rather than model 2 when we compare the parameters with the 
other PWNe, because model 2 of Kes 75 is clearly more extreme than model 1. However, 
we find Kes 75 peculiar in its results, even for the distance of 6kpc (model 1). The fitted 
values of ?7 = 5 X 10^^ and fpwN = 420km/sec are significantly smaller than other PWNe 
we studied so far. The parameters r] and wpwN become large, if we increase the local ISRF 
energy density f/isRF in the same manner as for model 2 of G54. 1+0.3. However, to get 
two orders of magnitude larger value of 77 and fpwN ~ lOOOkm/sec, more than an order 
magnitude larger Uisrf is needed (see Section 17.11 in details) and the age of the pulsar 
tage becomes around 300 yr, which is as small as the age of SNR Cassiopeia A. Because 
Kes 75 is the youngest PWN in our study, more precise studies of how a PWN inside a 
SNR is created may be important besides the magnetar-like properties of its central pulsar 
PSR J1846-0258. As for the particle injection, the parameters are not unusual except for a 
little larger value of 7b than other PWNe. Especially, the high energy power-law index at 
injection p2 = 2.5 is similar to other PWNe. 



Bucciantini et al. 



( I2OIOI ) studied Kes 75 with their spectral evolution model. There are 
two main differences between their model and ours. First, while we model the magnetic 
field evolution assuming energy conservation (Equation ([5])), they consider the adiabatic 
loss of the magnetic energy. Second, they consider dynamical evolution of a PWN inside a 
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SNR. The age of the system is determined from the dynamical properties of system in their 
model, while we determine it from the spectral properties. They assume that the age of 
the system is 650yr old, which is close to our result 700yr. The parameters of the particle 
injection pi = 1.7 and 7b = 8 x 10^ are almost consistent with our model pi = 1.6 and 
7b = 2 X 10^, but the high energy power-law index at injection p2 = 2.3 is a little harder 
than our value p2 = 2.5. The lower value of p2 increases the calculated flux in the hard 
X-rays, but their fitted spectrum still underpredicts hard X-ray observation. Their current 
magnetic field strength ~ 30yuG is almost consistent with our model -Bnow = 20yuG, but 
the adopted energy density of the ISRF ~ 24eV/cm^ is very large. However, because the 
assumed temperature of the ISRF Tisbf = lOOOK is high, Klein-Nishina effect significantly 
suppresses the IC/ISRF flux in their model. 



Bucciantini et al. 



torn . Only for 



The uniqueness of Kes 75 PWN also appears in 
the case of Kes 75, the fraction of the magnetic energy injection -Emag of the spin-down 
power L{t) (corresponding to rj in our model) is almost two orders of magnitude lower 
than other young PWNe in both studies. This relative smallness of t] compared with other 



young PWNe is like, 
is different between 



y to be real and i s very interesting, although the absolute value of rj 



Bucciantini et al. 



(120101 ) and the present work. The difference of the 



absolute value of rj is most probably due to the difference of the magnetic field evolution 
model. 



Lastly, we discuss the hard X-rays. 



McBride et al. 



(120071 ) argued that the hard X-ray 



emission detected by INTEGRAL is dominated by the emission from the PWN, i.e., it is not 
from the pulsar. However, it is difficult to reproduce the observed hard X-ray emission with 
the current magnetic field -Bnow = 20/iG. If bot h the observed soft an d hard X-ray spectra 



are fitted by power-law spectra as discussed in iMcBride et al 



cooling break frequency should be above lO^^Hz {B^ 



( 120071 ). the synchrotron 
5fiG for tagc = 0.7kyr) to 
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reproduce the observed soft X-ray photon index Tx < 2. This is difficult to reahze for the 
same reason as for G21.5-0.9 because the observed 7-rays to X-rays flux ratio constrains the 
current magnetic field strength. So that we think that the hard X-ray emission may have a 
pulsar origin. 



6. GO.9+0.1 



G O. 9-1-0.1 is a compo site super nova remnant and it s PWN componen t is observed in 



radio (IDubner et al. 



20081 ). X-rays (iGaensler et al. 



200lh. and TeV o'-ravs f|A 



2005 ). In X-rays, a jet-torus like structure is observed (IGaensler et al. 



2001 



laronian et al. 



Porquet et al. 



20031 ). In radio, GO.9+0.1 con sists of a compact 



outer SNR shell ^ 4' in radius ( 



Dubner et al 



shell is observed (iPorquet et al 



'WN which has a radius of ~ 1' and an 



20081 ) . No X-ray counterpart of the outer SNR 



20031). The central pulsar of GO.9+0.1 (PSR J1747-2809) 
has recently been detected in radio, with a period P = 5.22 x lO^^sec , its time derivative 



P = 1.56 X lO^^^sec ■ sec"^ (r^ = 5.3kyr) and unknown braking index (IGamilo et al.ll2009l ). 
So that we assume n = 3. The actual age of GO.9+0.1 is unknown. We consider two 
cases for the distan ce of GO.9+0.1: 8kpc (lower limit) and 13kpc (upper limit) given by 



Camilo et al. 



(120091 ) as model 1 and 2, respectively. The ISRF energy density of model 1 is 
larger than model 2, especially in optical band. We assume that GO.9+0.1 PWN is a sphere 
of radius 2.3pc at (r, z) ~ (Okpc, 15pc) or 3.8pc at (r, z) ~ (5kpc, 20pc) for model 1 and 2, 
respectively. The energy density of the ISRF is (t/m, Uqpt) = (1.6eV/cm^, 15eV/cm^) for 
model 1 and (?7ir, ?7opt) = (1.2eV/cm^, 2.0eV/cm^) for model 2, respectively. 



6.1. Model 1 



Figure [To] shows the model spectrum of GO. 9+0.1 when the distance to GO. 9+0.1 is 
taken to be 8kpc together with observational data. We fit the data with the parameters 
r] = 3.0x 10-3, = 2.0kyr, 7^,, = 8.0 x 10^ 7b = 7min = 4.0 x 10^ and p2 = 2.6. 
The observed radio spectrum is not fitted with the usual value of pi in a range of > 1 
but we take 7min = 7b- The fitted fraction parameter is close to that of the Crab Nebula, 
being just 0.6 times the value of the latter. The fitted age tage = 2.0kyr deviates from the 
characteristic age Tc = 5.3kyr, i.e., the current period is close to the initial period. The 
expansion velocity fpwN = 1120km/sec is comparable with that of G21.5-0.9 and model 
2 of G54. 1+0.3. The pulsar parameters are tq = 3.2kyr, Lq ■ tq = 1.2 x lO^^erg and the 
current total energy i?tot(2.0kyr) = 4.4 x lO^^erg, which is almost a half of Lq ■ tq. The 
current magnetic field strength of GO. 9+0.1 turns out to be Snow = 15fiG. 

A single power-law injection well describes the observations. The radio emission comes 
from the particles suffering from the adiabatic cooling and details will be discussed in 
Section [7711 The 7-ray emission is dominated by the IC/IR and IC/OPT. The synchrotron 
cooling break frequency z/c ~ 5 x lO^^Hz corresponds to the flux peak. 

Figure [11] shows the evolution of the emission spectrum (left panel) and the particle 
distribution (right panel) of GO. 9+0.1 with the use of the same parameters as in Figure [TTl 
While the synchrotron flux decreases, the inverse Compton flux increases with time. This 
feature is the same as G21.5-0.9 because a spin-down time tq = 3.2kyr is large. In the right 
panel, we see how a single power-law distribution evolves with time, producing the low 
energy tail. 
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6.2. Model 2 

Figure [12] shows the model spectrum of GO. 9+0.1 when the distance to GO. 9+0.1 is 
taken to be 13kpc together with observational data. We fit the data with the parameters 
r] = 1.0x 10"3, t,g, = 4.5kyr, 7^,, = 1.0 x 10^, 7^ = 7^;^ = 1-0 x 10^ and p2 = 2.6. The 
fraction parameter is similar to that of model 1. The main difference from model 1 is in 
the fitted age of GO. 9+0.1, which becomes more than twice that of model 1 because we 
need a larger current total energy Etot is needed given the larger distance to the object. 
The expansion velocity fpwN = 830km/sec is a little smaller than model 1. The pulsar 
parameters of GO. 9+0.1 are tq = O.Skyr and Lq ■ tq = 4.8 x lO^^erg. The current total 
energy ii^tot(4.5kyr) = 4.2 x 10^^ is an order of magnitude larger than model 1. Note that 
tage is larger than tq for model 2 and vice versa for model 1. The current magnetic field 
strength -Bnow = 12/iG is similar to model 1. 

The 7-ray emission is dominated by the IC/IR component. Although the energy 
density of the ISRF is very different from model 1, the energy density of the magnetic field 
remains similar, because of the strong Klein-Nishina effect on the IC/OPT in model 1. 
Note that the difference of the ISRF energy density between model 1 and 2 is mainly that 
in optical band. The synchrotron cooling break frequency is i^c ~ 2 x lO^^Hz. 

Figure [13] shows the evolution of the emission spectrum (left panel) and the particle 
distribution (right panel) of GO. 9+0.1 with the use of the same parameters as in Figure [13] 
(model 2). Both the synchrotron fiux and the IC fiux decrease with time, and this behavior 
is the same as G54. 1+0.3 shown in the left panel of Figure [5] This different evolution in the 
left panels of Figure [12] from in Figure [TT] arises from the difference between t^ge and tq in 
the models. 



- 26 - 



6.3. Discussion 

Both model 1 and 2 well reproduce the observed spectra and both models are 
acceptable. The models differ in the IC/OPT ffux, but the current observations cannot 
distinguish the models. A low energy, power-law component for the injected particles, is not 
needed (see the injection spectrum of the right panels of Figures [TT] and [T3l) . to reproduce 
the observed radio spectrum of GO. 9+0.1. This is because the adiabatic cooling of the 
injected particles can create the observed radio spectrum as will be discussed in Section [7^11 
However, we should note that no other PWNe which we have studied can be fitted by the 
single power-law injection of the particles. 



Fang fc Zhang] (|2010| ) studied the spectral evolution of GO. 9+0.1. They adopted the 
model in which GO. 9+0.1 is located at the galactic center {d = 8.5kpc in their model) 
with similar values of the ISRF energy densities f/jR = 0.5eV/cm'^ and Uopt = 20eV/cm'^ 
to our model 1. The most interesting point in their model is that the particle distribution 
at iniection is giv en by a relativistic Maxwellian plus single power-law distribution (e.g.. 



Spitkovsky 



20081 ). not a broken power-law distribution. They could reproduce the observed 
spectrum of GO. 9+0.1 with p2 = 2.5 and 7b ~ 4 x 10^. The current magnetic field strength 
8.1/iG is about a half of our value -Bnow = 15/iG. As will be discussed in Section m| their 
conclusion is consistent with our study, but the existence of the relativistic Maxwellian 
component is not essential to reproduce the observed radio spectrum. We think that 
adiabatic cooling of the injected particles could also works in their model to reproduce the 
observed radio spectrum. 
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7. DISCUSSIONS AND CONCLUSIONS 

We discuss first the dependence of emission spectrum on tlie adopted ISRF energy- 
densities. Next, we discuss about tlie difference and similarity of the fitted and derived 
parameters among five young TeV PWNe including the Crab Nebula. We search for 
correlations between the central pulsar properties and the fitted parameters. 



7.1. Discussion 

There is a possibility that the ISRF energy density is locally different from the mean 
values of the Galaxy as discussed in the case of model 2 of 054.1+0.3 and we see the 
7-ray emission is dominated by the IC/ISRF except for the Crab Nebula. To reproduce 
the observed power of the IC/ISRF, a larger ISRF energy density f/isRF leads to a smaller 
current total energy of particles (1 — i])Etot{ta,ge)- Accordingly, the fraction parameter 
7] needs to be larger to reproduce the observe d power of the synchrotron radiation. We 
showed it in our previous paper (Section 3.3 of iTanaka &: Takaharall2010( ) that the power of 



the synchrotron radiation and that of the IC/ISRF roughly behave as Pgyn oc (1 — ri)'r]E'^ 



tot 



and -Pic/isRF oc (1 — i])EtotUisRF, respectively. These relations lead to i?tot ^isrf 

1] oc f/ilj^F together with observed Psyn and Pic/isrf- As for tage, we find that a larger f/isRF 

leads to a smaller tage from the integration of Equation ([T]). 

However, Klein-Nishina effect makes the dependence of E^ot and r] on f/jsRF somewhat 
milder. As seen from Tabled the adopted f/iR in model 2 of G54.1+0.3 is four times larger 
than that in model 1, but Etot of model 1 is about three times larger than that of model 2. 



Accordingly, i] of model 2 is about seven times larger than that of model 1. In section [573| 
we considered how large f/isRF is required in order for t] of Kes 75 to be as large as other 
PWNe. We can estimate that more than an order of magnitude larger ?7isrf is required for 
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model 1 and almost two order of magnitude larger f/isRp is required for model 2. 



For the fraction parameter rj, all the young TeV PWNe have a value much smaller than 
unity and most of them are similar to each other. We conclude that the fraction parameter 
7] in our model is typically a few x 10^^ for young TeV PWNe except for Kes 75. 

For the break energy, 7b ~ 10^~^ is found for all the young TeV PWNe. Together 
with the minimum energy 7inin and the low energy power-law index at injection pi, these 
parameters determine the supply rate of the particles from the pulsar wind and thus 
determine the pair multiplicit y k inside the pulsar magnetosphere and the bulk Lorentz 



factor of the pulsar wind Pw (ITanaka fc Takahara 



2010 



Bucciantini et al. 



2010|). For the 



typical values of the power-law indices at injection 1 < pi < 2 < p2 < 3, the particle number 
flux is estimated as Nmj(t) ~ 7min<5inj(7min, ^) and the particle energy flux is estimated 
as Epartit) ~ 7b"^eC^<5inj(7b5 We can estimate k as A^inj = uNqj, where Nqj is the 
Goldreich- Julian number flux and Pw is estimated as L{t) ~ ii^part(^) ~ Ni^j{t)T^meC^. The 
fitted minimum energy (7min < a few xlO^) gives a lower limit of k and an upper limit of 
Pw However, the values of k and P„ of GO. 9+0.1 is fixed because 7min = 7b- Derived k and 
P„ for each PWN are listed in Table |2l We find that the lower limit of k, is larger than 10"^ 
for all young TeV PWNe and that the upper limit of P^ is smaller than 7b for the Crab 
Nebula and Kes 75 (model 1 and 2). These quantities may be more constrained by the 
future lowest frequency radio observations, such as LOFAR, ASKAP and SKA. 

The high energy power-law index at inje ction ^ 2.5 takes a v ery similar value 



Bucciantini et al 



( 120101 ). spatially variable 



for all five young TeV PWNe. As stressed by 
fluxes and photon indices of X-ray observations are not guaranteed to be reproducible 
in a one-zone broad band emission models. For example, the hard X-ray observation of 
Kes 75 may suggest that the hard X-ray emission mainly comes from the central region 
of strong magnetic fields. However, the result of p2 ~ 2.5 for all the young TeV PWNe 
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is still interesting. It suggests that acceleration process at the pulsar wind termination 
shock and/or the cool ing, advectioii and diffusion of the accelerated particles are common 



to young PWNe (e.g. 



Bambaetal 



20101 ) 



The low energy power-law index at injection pi, on the other hand, is different for each 
PWN and varies 1.0 < pi < 1.6. Moreover, in the case of GO. 9+0.1, we either do not need 
a low energy component or it should be very hard {pi ^ 1.0). On the assumption of an 
uniform PWN, this behavior of GO. 9+0.1 is explained in the following way. Low energy 
particles are mainly cooled by adiabatic cooling. When we take into account only the 
adiabatic cooling. Green's function of the continuity equation of particles in energy space 

-G(7, t) - a— (|g(7, t)) = 5(7 - loM - to) (8) 

becomes 

G(7, t) = - ■ 7~'0(7o - 7)5 f In f + - 1^ -) ' (9) 
a \ to a 7o/ 

where a takes into account an accelerated (a > 1) or decelerated (a < 1) expansion 
of the PWN (-RpwN oc t") and Q{x) is the Heaviside step function. For simplicity, we 
consider the solution of the time-dependent injection q{j,t) = qQ{t/ti)~^5{j — 7i)0(t — ti), 
where go = const, is a normalization factor, ti is the time when the injection starts, 71 is 
the particle energy at injection and /3 is the time dependence of injection. The particle 
distribution N{'y,t) is given by 

We consider the PWN expanding at a constant velocity a = 1 and the particles whose 
energy is lower than 71 = 7b = 7min for GO. 9+0.1. Equation (|TOl) gives the particle 
distribution oc 7"^ for 7 < 71 = 7min- From Equation ([T]), for t <^ tq, /3 is while /3 
is {n + l)/{n — 1) for t ^ Tq. However, for t^gc ~ To ~ O(kyr) as in the both models 
of GO. 9+0.1, (3 smoothly varies from to 2 for = 3 on the time scale of tq. As seen 
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in the right panels of Figures [TT] and [T3] (see thick dotted curves at 10 kyr), the particle 
distribution for 7 < 7b continuously changes fr om to —2 on the time scale of tq. Thus, 



20081 ). corresponding to the 



the observed radio spectral index = —0.18 (iDubner et al.l 
power-law index ~ —1.4 of particle distribution, can be almost reproduced by the adiabatic 
cooling. Note that the slope of the particle distribution for 7 < 7b is slightly different 
between the right panels of Figures [TT] and [T31 because Tq is different. Note also that the 
observed radio spectral index ar = —0.18 is clearly smaller than = 1/3 which corresponds 
to the low frequency tail of synchrotron radiation. 

Concerning the properties of the central pulsars, when the braking index n is given, the 
spin-down evolution of the pulsar is characterized by two parameters, the initial spin-down 
power Lq and the spin-down time Tq in Equation ([1]). In other words, the individualities of 
each pulsar come from these two quantities which theoretically represent the two parameters 
equivalent to the initial rotational energy and the magnetic energy of the pulsar. In Figure 
[m we plot the correlation between the initial rotational energy Lq ■ tq versus the magnetic 
energy of the pulsar Eb = B^Rl/Q, where oc pi/2pV2 (assuming magnetic dipole 
radiation) and i?* = lO^cm are the surface dipole magnetic field and the radius of the 
pulsar, respectively. Figure [H] shows an anticorrelation between Lq-tq and Eb, although the 
statistics is rather poor and Kes 75 dominates in this anticorrelation. We should mention 
that this anticorrelation remains unchanged even if we use the canonical value n = 3 for 
Kes 75, instead of we have used n = 2.65. Although the canonical value n = 3 keeps Eb 
unchanged from its birth, when n is not 3, it seems better to use the initial value of Eb 
(Ebq) rather than that at present time -EB(^agc)- However, for Kes 75 {n = 2.65), -EB(^agc) is 
only a factor of 1.3 smaller than Ebq- The value of n also changes Lq ■ tq. However, Lq ■ tq 
for n = 2.65 is only a factor of 1.4 larger than that for n = 3 with fixed -Etot(^age), which we 
can almost independently determine from the observed power of the IC/ISRF. 
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We also search for correlations between the parameters of the injection spectrum and 
the pulsar parameters. The pulsar parameters include the spin-down power L{t) oc P/P^, 
the surface dipole magnetic field oc P^/^P^/^, the potential difference at the polar cap 
<l>poie oc pyyp^/^ and the light cylinder magnetic field Pic oc B^{R^VL/cf oc pi/2/p5/2^ 
where VL is the current angular velocity of the pulsar. In the left panel of Figure [151 we plot 
the correlation between the fraction parameter r] versus the light cylinder magnetic field 
Pic. Although it seems to show some correlation, when we ignore Kes 75, the correlation 
is insignificant and the values of rj for the other objects spread only in a range 10~^ - 
10~^. In the right panel of Figure [TSl we plot the correlation between the maximum 
energy 7max versus the potential difference at the polar cap $poie- Because e'l'poie gives the 
maximum available electri c energy of the pulsar to accelerate particle, this correlation is 



expected as mentioned by 



Bucciantini et al. 



(120101 ). However, we do not find a significant 



correlation. This may be partly because only an upper limit of 7max is obtained, except for 
the Crab Nebula. For other combinations of the parameters, we do not find any significant 
correlations and we do not show them here. 

Lastly, we discuss about the age tage and the expansion velocity wpwN of PWNe. Our 
spectral evolution model can estimate the age of the central pulsar in a fairly reliable way 
from the observed power of the IC/ISRF. In contrast, the characteristic age Tc of a pulsar 
is not guaranteed to match the age of young pulsars, when the spin-down time tq is close 
to the age of the pulsar tage as seen in Equation ([2]). For example, Tc of the central pulsars 
inside the Crab Nebula and G21.5-0.9 are 1.2kyr and 4.8kyr, but tage fitted with our model 
and the ages estimated from the observed expansion rate of the synchrotron nebula are 
both ~ l.Okyr. 

We discuss about the properties of the SN explosion which creates the pulsars with 
obtained fpwN and Ptot- Simply, we expect the relation i'pwn ~ (-E'tot(^)/-E'sN)°'^VsN, where 
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EsN is the energy of SN ejecta and Vsn is the velocity of the front of the freely expanding 



ejecta. We derive this relation from the radius of PWIN 
ejecta estimated as -Rpwn(^) ~ (-E'tot(^)/-SsN)'''^VsN^ by 



e -RpwN in a freely ex p andin g SN 



van der Swaluw et al 



teoOlh . We 



do not find the relation fpwN oc -Etot > which suggests that the values of E^n and Vsn are not 



common to the PWNe. On the other hand, because the derived values of fpwN and E\ 



0.2 

tot 



differ by a factor of a few, the combination of Vsn / E^^ also ranges within a factor of a few. 



7.2. Conclusions 

In this paper, we apply our spectral evolution model to four young TeV PWNe, 
G21.5-0.9, 054.1+0.3, Kes 75, and GO. 9+0.1. We have succeeded in reproducing many 
observed properties of these PWNe based on this rather simplified one-zone model. 

The current observed spectra of all four young TeV PWNe are reconstructed with 
small values of the fraction parameter r/ ^ 1 as well as the Crab Nebula, i.e., the magnetic 
energy of these PWNe accounts for a very small fraction of the current total energy injected 
into the PWN -Etot(^age)- The fitted fraction parameters are typically a few xlO~^ except 
for Kes 75. The fraction parameter of the peculiar object Kes 75 is more than two orders of 
magnitude smaller than the typical value. 

The TeV 7-ray emission from the young TeV PWNe is dominated by IC/ISRF. Since 
the energy density of the local ISRF around the objects is somewhat uncertain, it is 
important to take into account its effect as considered in model 2 of G54. 1+0.3. On the 
other hand, the 7-ray emission at early phase of their evolution (e.g., tage < 300yr) is always 
SSC dominant because the magnetic energy density of the PWN is much larger than the 
local ISRF energy density. 

A broken power-law injection of particles well reproduces the observed spectrum from 
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radio through TeV 7-ray, except the case of GO. 9+0.1 where we do not need the low energy 
component. The observed spectrum of GO. 9+0.1 in radio is created by the adiabatic coohng 
of the high energy component of the injected particles. 

The fitted break energy 7b ~ 10^~^ is rather common. On the other hand, the fitted 
maximum energy 7max is a lower limit (more than ~ 10^) and the fitted minimum energy 
7min is an upper limit (less than ~ 10^). The high energy power-law index at injection 
P2 ~ 2.5 is common for all young TeV PWNe, while the low energy power-law index at 
injection pi varies in the range 1.0 < pi < 1.6. 

The fitted parameters of the injection spectrum in Equation gives a lower limit of 
the pair multiplicity n, which turns out to be more than 10^. We have also estimated an 
upper limit of the bulk Lorentz factor of the pulsar wind r„. For G21.5-0.9 and G54. 1+0.3, 
the upper limit of is still consistent with 7b, but the obtained F„ of Kes 75 is clearly 
less than 7b. The latter feature is the same as for the Crab Nebula in our previous work 
Fw < 10^ < 7b ~ 10^. On the other hand, for GO. 9+0.1, the value of F„ is similar to 7b 
without uncertainty in our model. 

From the fitted age, we can derive the initial spin-down luminosity and the spin-down 
time of the central pulsar. We take the initial rotational energy Lq ■ tq and the magnetic 
energy E-q of the pulsar as two independent parameters which characterize the spin-down 
evolution of the pulsar and search for a correlation between them. They seem to be 
anticorrelated, although the statistics is rather poor. We also search for correlations 
between the fraction parameter and the parameters of the injection spectrum versus the 
central pulsar properties. However, we find no significant correlations. 

We are grateful to Y. Ohira for useful discussions. We thank the anonymous referee for 
his/her meticulous reading of the manuscript and helpful comments. This work is partly 
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Fig. 1. — Model spectrum of G21.5-0.9 at tagc = l-Okyr for model 1, where fitting to the 
infrared observation is included. The solid line is the total spectrum which is the sum of the 
synchrotron (thin-dashed line), IC/CMB (dotted line), IC/IR (dot-dashed line), IC/OPT 
(dot-dot-d ashed line) and SS C (dashe d line) spectra, r espec tively. The observed data ar e 
taken fromiSa. 
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Fig. 2. — Evolution of the emission spectrum (left panel) and the particle distribution (right 
panel) of G21.5-0.9 for model 1. The solid line, dashed line, dot-dashed line, and dotted line 
correspond to 300yr, Ikyr, 3kyr, and lOkyr from birth, respectively. The thin dotted line in 
the particle distribution (right panel) is the total injected particles at an age of lOkyr. 
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Fig. 3. — Model spectr um of G21.5-0.9 at tp,„ p, = l.Okyr for model 2, where we ignore 
the infrared observation (iGallant fc Tuffslll998l ) in the spectral fitting. Used and obtained 
parameters are tabulated in Table [1] 
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Fig. 4. — Model spectrum of G54.1+0.3 at tago = 2.3kyr for model 1. The observational data 
and the lyr, 5cr sensitivity for Fermi LAT are plotted. This sensitivity is very similar to the 
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Fig. 5. — Evolution of the emission spectrum (left panel) and the particle distribution (rij 
panel) of G54. 1+0.3 for model 1. 
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Fig. 6. — Model spectrum of 054.1+0.3 at tago = l-7kyr for model 2, where an enhanced 
local ISRF is assumed. Fermi LAT sensitivity and the observed data are the same as in 
Figure HI Used and obtained parameters are tabulated in Table [1] 
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Fig. 7. — Model spectrum of Kes 75 at 



age 



v[Hz] 

= 0.7kyr for model 1, where the distance is 



plotted. The observe d data are taken from 
(radio), 



Helfand et al 



020 031); 



Salter et al. 



Morton et al.i fi2007i ) (X-ray), 



(1989 


); 


Bock & Gaens 


lei 


( 


2005) 


. Djannati-Atai' et al. 


(200" 


n h- 



ray). Used and obtained parameters are tabulated in Tabled! 



-47- 




Fig. 8. — Evolution of the emission spectrum (left panel) and the particle distribution (rij 
panel) of Kes 75 for model 1. 
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Fig. 9. — Model spectrum of Kes 75 at tage = 0.88kyr for model 2, where the distance is 
taken to be lO.Gkpc. Fermi LAT sensitivity and the observed data are the same as in Figure 
[3 Used and obtained parameters are tabulated in Table [1] 
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Fig. 10. — Model spectrum of GO.9+0.1 at tago = 2.0kyr for model 1, where the distance 
is taken to be 8kpc. The observational data and the lyr, So" sensitivity for Fermi LAT are 
plotted. However, this sensitivity curve may be more worse because GO.9+0.1 is significantly 
closer to the Galact ic center than G21.5-0.9 on th e sky. T he observed data are ta ken from 
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Fig. 11. — Evolution of the emission spectrum (left panel) and the particle distribution 
(right panel) of GO. 9+0.1 for model 1. The total injected particles at an age of lOkyr (thin 
dotted line in right panel) shows that the injection spectrum is given by a single power-law 
distribution. Note that the pileup feature is appeared in the right panel of the particle 
Lorentz factor 7 < 10^ for t — 0.3kyr. This is made of numerical error and the particles in 
this feature do not contribute the emission spectrum oi u > lO^Hz. 
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Fig. 12. — Model spectrum of GO. 9+0.1 at tagc = 4.5kyr for model 2, where the distance is 
taken to be 13kpc. Fermi LAT sensitivity and the observed data are the same as in Figure 
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Fig. 13. — Evolution of the emission spectrum (left panel) and the particle distribution (right 
panel) of GO. 9+0.1 at 13kpc. As seen from the right panel, the injection spectrum is a single 
power-law distribution. Note that the pileup feature is also appeared in the right panel of 
the particle Lorentz factor 7 < 10^ for t = O.Skyr and is the same as the right panel of Figure 
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Fig. 14. — The correlation between the initial rotational energy Lq ■ tq versus the magnetic 
energy B?^^R^^^^^/6 of the central pulsars. Model 1 of G54.1+0.3 and model 2 of Kes 75 are 
not plotted. 
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Fig. 15. — The correlations between the fraction parameter r) versus the magnetic field 
strength at the light cylinder Bi^ (left panel) and the maximum energy 7max (lower limit 
except for the Crab Nebula) versus the potential difference at polar cap $poie (right panel). 
Model 1 of G54. 1+0.3 and model 2 of Kes 75 are not plotted. 



Table 1. Adopted parameters and fitted and derived parameters to reproduce the current observed spectrum. 



Symbol 


Crab" 


G21.5-0.9 


G21. 5-0.9*' 


G54.1+0.3 


054.1+0.3'^ 


Kos 75 


Kos 75'* 


GO. 9+0.1 


GO.g+O.l"* 


Model 




1 


2 


1 


2 


1 


2 


1 


2 



Adopted Parameters 



d{kpc) 


2.0 


4.8 


4.8 


6.2 


6.2 


6.0 


10.6 


8.0 


13 


-RpWN,now(P'^) 


1.8 


1.0 


1.0 


1.8 


1.8 


0.29 


0.50 


2.3 


3.8 


P(msGc) 


33.1 


61.9 


61.9 


136 


136 


326 


326 


52.2 


52.2 


P(10~^^) 


4.21 


2.02 


2.02 


7.51 


7.51 


70.8 


70.8 


1.56 


1.56 


Ti 


2.51 


3.0 


3.0 


3.0 


3.0 


2.65 


2.65 


3.0 


3.0 


(7lR{cV/cm=») 




1.0 


1.0 


0.5 


2.0 


1.2 


1.0 


1.6 


1.2 


(7oPT{oV/cm^) 




2.0 


2.0 


0.5 


0.5 


2.0 


2.0 


15 


2.0 


Fitted Parameters 


r,(10"3) 


5.0 


15 


8.0 


0.3 


2.0 


0.05 


0.006 


3.0 


1.0 


tago (kyr) 


0.95 


1.0 


1.0 


2.3 


1.7 


0.7 


0.88 


2.0 


4.5 


7max(103) 


7.0 


> 2.0 


> 2.0 


> 1.0 


> 1.0 


> 1.0 


> 0.8 


> 0.8 


> 1.0 


7b(10^) 


6.0 


1.2 


0.7 


3.0 


1.8 


20 


50 


0.4 


1.0 


7mi„{10^) 


< 0.1 


< 3.0 


< 3.0 


< 20 


< 20 


< 5.0 


< 5.0 






Pi 


1.5 


1.0 


1.0 


1.2 


1.2 


1.6 


1.4 






P2 


2.5 


2.55 


2.5 


2.55 


2.55 


2.5 


2.5 


2.6 


2.6 


Derived Parameters 




1800 


980 


980 


770 


1040 


420 


560 


1120 


830 


S„ow{mG) 


85 


64 


47 


6.7 


10 


20 


24 


15 


12 


To(kyr) 


0.7 


3.9 


3.9 


0.6 


1.2 


0.2 


0.003 


3.2 


0.8 


Lo . ToClO^^^erg) 


74 


6.5 


6.5 


5.4 


2.6 


1.5 


210 


12 


48 



^Results are taken from lXanaka Sz Takaharal J201[j'l ■ 

'^All the adopted parameters are the same as model 1 of G21.5-0.9, but ignoring the observation in infrared. 
'^Assumed Uij^i is different with model 1 of G54. 1 + 0.3. 

■^Assumed distances to the objects are different with model 1 of Kes 75 and GO. 9+0.1, respectively. 



Table 2. The derived pair multiplicity and bulk Lorentz factor. 



Symbol 


Crab 


G21.5-0.9 


G21.5-0.9 


G54.1+0.3 


G54.1+0.3 


Kes 75 


Kes 75 


GO.9+0.1 


GO.9+0.1 


Model 




1 


2 


1 


2 


1 


2 


1 


2 


k(10*) 


> 420 


> 13 


> 19 


> 3.7 


> 5.2 


> 2.8 


> 0.85 


8.3 


3.4 


rw(io5) 


< 0.07 


< 0.67 


< 0.53 


< 2.1 


< 0.91 


< 1.4 


< 4.6 


1.1 


2.7 



